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Abstract 

The ageing Lie algebra oge(d) and especially its local representations for a dynamical 
exponent z = 2 has played an important role in the description of systems undergo- 
ing simple ageing, after a quench from a disordered state to the low-temperature phase. 
Here, the construction of representations of age(d) for generic values of z is described 
for any space dimension d > 1, generalising upon earlier results for d = 1. The mecha- 
nism for the closure of the Lie algebra is explained. The Lie algebra generators contain 
higher-order differential operators or the Riesz fractional derivative. Co-variant two-time 
response functions are derived. Some simple applications to exactly solvable models of 
phase separation or interface growth with conserved dynamics are discussed. 



1 Introduction 



Understanding the cooperative properties of strongly interacting many-body problems continues 
to pose many challenges. Here, we are interested in a particular class of non-equilibrium phe- 
nomena, usually referred to as ' ageing'. While first systematically studied in structural glasses 
quenched from a molten state to below the 'glass-transition temperature' [51], very similar 
phenomena have been found in many other glassy and non-glassy systems far from equilib- 
rium, see e.g. [H |22] for surveys and further references. Schematically, one may characterise 
ageing systems by (i) a slow relaxation dynamics, (ii) absence of time-translation- invariance 
and (iii) dynamical scaling. Therefore, if ageing is understood in this way, the study of the 
ageing of systems far from equilibrium gives physically well-motivated setting for the analysis 
of the consequences of dynamical scaling and the investigation of possible generalisations and 
extensions. 

Throughout this paper, we shall restrict to what is known as 'simple ageing', where the 
dynamics can described in terms of a single time-dependent length L(t) such that for large times 
Lit) ~ t l l z , which defines the dynamical exponent z. The analysis of the ageing of several simple 
magnetic systems, without disorder nor frustrations, without any macroscopic conservation law 
of the dynamics, and undergoing ageing when quenched to a temperature T < T c below the 
critical temperature T c > is characterised by the dynamical exponent z = 2 [6] . Furthermore, 
it has been suggested that the detailed scaling form of the two-time correlators and responses 
can be obtained by an extension of simple dynamical scaling with z = 2 towards a larger Lie 
group PJ2, [16j [T7] of which the so-called ' ageing algebra' agt(d) = (X i, Y$. , M , Rij) 

2 / \<i<j<d 



is the Lie algebra. It can be defined in terms of its non-vanishing commutators 

[Xn,Y^] = Q- m )y» m , [X n ,X n ,] = (n-n')X n+n ,, \Y® ,Y^] = ^M , 

[Rij, Rke] = 5 ie R jk + 5j k Ru — 6 ik Rj£ — djeRik, [Rij, Yj£'] = d~ jk Y^ — 5 ik Y$ (l.l) 

with n, n' = 0,l,m = ±i and 1 < i < j ' < d. When acting on space-time coordinates (t, r), a 
representation of (II. ip in terms of affine differential operators is: 

x 

X = —td t — -(r ■ d r ) — — dilatation 

X\ = —t 2 d t — t(r ■ d r ) — — r 2 — (x + £)t special transformation 
Y_\ = —d n space-translations (1.2) 

2 

Y^ = —td n — Airi Galilei-transformation 

2 

M Q = —M phase shift 

Rij = Tid Ti — Tjd ri = —Rji, rotations. 

together with the physical interpretation of these generators. Clearly, the representation fl 1.2ft 
has a dynamical exponent z = 2 and acts locally on the space-time coordinates. Furthermore, 
it generates a set of dynamical symmetries of the Schrodinger (or diffusion) equation: 

S<f){t, r) = \ 2Md t + ^{x + £ - d/2) - V^j <j>(t t r) = 0, (1.3) 
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in the sense that each of the generators of Oflc(d) maps a solution of (II. 3p onto another solution. 
Herein, (Ai,x,^) is a triplet of constants which together characterise the solution = (f)(M,x,£) 
of this equation. Physically, Ai is interpreted as an inverse diffusion constant, whenever A4. G M. 
or else as a non-relativistic mass if M. = im with m G R. Furthermore, x and £ are the two 
independent scaling dimensions which characterise the field = 0(.m,z,5)lJ 

For systems undergoing simple ageing with z — 2, the extended local scale-invariance as 
described by the representation ( II. 2p of agc(d) indeed seems to give an appropriate description, 
including several exactly solved examples where £ ^ is required [332 [22] - the best-known 
example is the ID Glauber-Ising model quenched to T = 0. The main prediction concerns 
the derivation of two-time (linear) responses R = R(t, s) = S gf$) of the order parameter 

with respect to its conjugate magnetic field, from the assumption of co- variance under the 
chosen representation of agc(d). The derivation is based on the Janssen-de Dominicis theory 

[27] which permits to re-express responses R(t,s) = ^0(t)0(s)^ as a correlator of the order 

parameter = (fi(M,x£) and the conjugate 'response field' = 4>t M g|S (see e.g. [27J, [531 E3] for 
introductions and detailed references). Co- variance under a$je(<f) then leads to a set of linear 
partial differential equations for R(t, s)H In statistical physics, a common formulation uses a 
stochastic Langevin equation 

d t( j>(t,r) = -D^L + v (t,r) (1.4) 

with a Ginzburg-Landau functional "H and a centred gaussian noise r\ with a 5-correlated second 
moment. In the standard Janssen-de Dominicis formalism, this is related to the equation of 
motion derived from the dynamic functional 



J[0,0]= dtdr 



\2Md t - A r ) + ^jP- - T<f) 2 



l - Jdrdr'Mr)a(r-r')Mr') (1-5) 



where D 1 = 2.M, V is the 'potential' part of the Ginzburg-Landau functional, a(r) the 
correlator in the initial state 4>o(r) = 0(0, r) and T the temperature of the external heat bath. 
In the case of a quadratic potential V, this action produces two equations of motion which 
both have a form similar to (11.21) . and with the above characterisations of and [42]. In this 
case, one has a decomposition ^"[0,0] = ^[0,0] + Jb[4>} such that the 'deterministic part' Jq 



1 If one were to add the time-translations X-i = ~dt, one would obtain the embedding agc(d) C age(cZ) © 
CX-i sct)(d) into the Schrodinger algebra scf)(d), with (jl.ljl extended to n = —1. In the representation (|1.2I) 
one then has £ = 0. Invariance of (|1.3|) implies further that a; = d/2, a fact essentially known to Jacobi and Lie 
in 1843 and 1881, respectively and since then re-discovered many times. 

Since time-translation-invariance is absent by definition in ageing systems and hence a stationary state is 
never reached, the new scaling dimension £ is a further universal characteristics of the relaxation process. The 
transformation <fi(t, r) = i^*<i>(t, r) maps an age(c?)-quasi-primary scaling operator <j> with the scaling dimensions 
{ x <P>£,<t>) to the 5cf)(d)-quasi-primary scaling operator $ with scaling dimensions = + 2^,^$ = 0) and 
confirmed for two-point functions |19| . The relationship between three-point and higher correlators is more 
involved [38] . 

2 There is a systematic way to 'dualise' the parameter M. which permits to further extend ctge(d)-covariance 
to the covariance under a parabolic sub-algebra of the conformal algebra conf(c? +2) in d + 2 dimensions, 
which gives an algebraic way to derive the causality condition t > s necessary for the interpretation of the co- 
variant two-point function as a physical response [HI [H] • In this dualised form, the non-relativistic AdS / CFT 
correspondence becomes evident, see e.g. [Tl 1571 14*5] . 
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is invariant under the action of the Galilei sub-algebra gal(d) = \Y_fl, M , Rij) . This 

\ ± 2 / l<i<j<d 

implies the Bargman super-selection rules [2] . From it follows that pT2j [22] 

1. the computation of all response functions can be exactly reduced to the computation of 
response in the 'noiseless' theory governed by the functional j7o, whose averages are given 
by co-variance conditions under the dynamical symmetry algebra 

2. all correlators can be exactly reduced to certain integrals over higher response functions 
of the 'noiseless' theory 

These calculations have been carried out for a long list of models undergoing simple ageing 
with z — 2 p$J [HI [EH [22]. In spite of many encouraging numerical results, there is not yet 
a satisfactory analytical treatment available for a non-quadratic potential V, where the main 
difficulty come from the non-trivial co-factors of the scaling operators <fr, <fi under Galilei- and 
special transformations (Yi, X±). A way around this problem might be to perform first a 
dualisation with respect to Ai, which converts the projective representation (jl.2p of the non- 
semi-simple Lie algebra agc(d) into a true representation of a sub-algebra of conf(c/ + 2) [19| I5U]. 

Can one extend this procedure, at least for linear stochastic Langevin equations of motion, 
to arbitrary values of the dynamical exponent z ? If we were to restrict to locally realised 
algebras, the recent classification of the non-relativistic limits of the conformal algebra [TD] 
would only admit the cases (i) z — 1: the conformal algebra cotlf(d) or the conformal Galilean 
algebra CGA(d) [HI EH SO] , eventually with the exotic central extension for d = 2 [32] (ii) z = 2: 
the Schrodinger algebra and (hi) z = oo; along with their sub-algebras. Further examples can 
only be found when looking at non-local representation, which means that we must go beyond 
the setting for first-order linear (afline) differential operators in the generators and must look 
for different, non-local realisations of the known abstract algebras. Some partial information is 
already available to serve as a guide: 

1. the Galilei-invariance of the non-relativistic equation of motion S<p = should be kept, 
which means that algebraically one should require 

[if \ Y°?] = %M , [S, Y^i] = \ { £S (1.6) 

2 2 - S1 2 

where the scalars A± are to be determined. For example, in sub-representations where 
Mq = 0, this can be extended to a Schrodinger algebra representation for any value of z 
[3D] , which describes massless particles, but with a rather trivial equation of motion, see 
below. 

2. In the context of ' local scale-invariance' different realisations of generalised symmetry 
algebras have been constructed by using certain fractional derivatives [TTJ [2U [22] • The 
closure of these sets of generators can only be achieved by taking a quotient with respect 
to a certain set of 'physical' states. Although this has been successfully applied to certain 
physical models [U [9] the closing procedure is not completely determined and it is not 
clear how to obtain the group (finite) transformations. 

3. There are many physically well-motivated systems, such as the Kardar-Parisi- Zhang equa- 
tion in ID with a dynamical exponent z = 3/2 [22] or else diffusion-limited erosion [21] 
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with z = 1 and a non-local equation of motion. One would like to be able to discuss their 
dynamical symmetries in a systematic way, but these cases do not seem to be included 
in the list of known representations of known algebras. 

A distinct and potentially more promising method has been explored in [23]. Therein, new 
non-local representations of age(l) for an integer- valued dynamical exponent z = n £ N were 
constructed. They read 



r 2 - \v* 



(1.7) 

The commutation relations ( II. ip are indeed satisfied, but with a notable exception, namely 

[Xx.yj] = gear's, (i.8) 

where the 'Schrodinger equation' takes the form 

SiP(t, r) := (zfid t - d z r )iP(t, r) = 0, (1.9) 

Therefore, the algebra is now closed only in the quotient space over solutions of eq. (II. 9p and 
in addition the representations (II. 7p act as dynamical symmetries of the same equation [23] . 



x = 


tO t — 

2 


1 a x 
-ro r 

2 2 


x x = 


- 


- trd r - (x + g)t) 


Y i = 

2 


— d r 




Yi = 


-td?- 1 - 


fir 


Mo = 


-li. 





In this paper, we shall show how to generalise this construction (II. 7p to any spatial dimension 
d > 1, which will be important in order to apply it to many physically relevant models. In 
contrast to the local case z — 2, the passage to d > 1 is not trivial. We shall describe this in 
section 2. Co- variant two-point functions are computed from these non-local representations in 
section 3, first for the case when z is even (with special attention given to the case z = 4) and 
afterwards for generic values of z, when the Riesz fractional derivative will play an important 
role. In section 4, we shall apply these results to some simple physical models, namely the 
kinetic spherical model with a conserved order-parameter and quenched to T = T c and the 
Mullins-Herring (or Wolf- Villain) equations of interface growth with mass conservation. The 
full space-time response are calculated from the non-local representations of agc(d) and the 
results will be compared with the known exact results [301 EH S3 H] • We conclude in section 5. 



2 Non-local representations of a$c(d) 

As we shall now show, the non-local representations of agc(d) are quite different for even and 
generic values of the dynamical exponent z. We shall discuss these cases separately. 



2.1 Representations of agt(d) with an even dynamical exponent 

Consider the following representations of agt(d) for an even dynamical exponent z = In: 
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Y {i) ■= -d 

2 -1/2 ' Ul 

Y% := -tdtA"- 1 - fin 

X 

X := -ntd t - -(r ■ d r ) - - 



V 2 
—r 

2 

M := -n {2.1) 



X 1 := {-nt 2 d t -t{r-d r )-{x + Ot)^ l -^r 2 



where di = d/dri, i = 1,2, ... ,d and A = A r = <9j<9j is the Laplacian. The spatial rotations 
Rij are unchanged with respect to ( II. 2p . In order to check the commutation relations ( II. ip we 
use the following relations: 

[A n ,rJ = 2nd i A n ~ 1 

[A n , r 2 } = Anr ■ d r A n ~ l + 2n(d + 2(n - 1)) A n_1 . (2.2) 
Then the commutators (II. ip are indeed satisfied, except the following: 

[X l ,Y®] = {n-l)1?d i A n - 2 S, (2.3) 

2 

where 

5 = 2nM-A" + f (x + f -^±^), (2.4) 

that is the generators ( 12. ID closed into the Lie algebra (XQc(d) only over functions from the 
solution space of the "Schrodinger" equation 



S*(t, r) = (2nfid t - A" + ^ fx + £ - - + ^ - \ \ $(t, r) = 0. 



(2.5) 



The representation (12.11) acts as a dynamical symmetry of the same equation. Since (I2.5P is 
linear, it is enough to check the following commutators 

[S, Y {t l] = [S, Y®] = [S, M ] = [S, R l3 ] = 

2 2 

[S,X ] = -nS, [S,Xx] = -2tnA n " 1 5. (2.6) 

Hence for any solution of 5$ = 0, the infinitesimal change 5$ = XS§ with a generator X 
chosen from (12. ip also solves the same equation. 



2.2 Representations of agc(d) with a generic dynamical exponent 

For generic values of the dynamical exponent z, an explicit representation of agc(d) is more 
difficult to find and turns out to be quite distinct from those found for an even- valued dynamical 
exponent. In particular when the dynamical exponent is an odd number a generalisation 
of one dimensional case in manner similar to z-even case is possible only for a sub-algebra 
{X , Xx, Y\, Mq\ which does not include the rotations. In addition this algebra is closed over 
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solution of the system of d equations, which are not invariant under entire algebra. We shall 
not develop this case here. 

Instead, in order to treat the case with a generic value of z, recall first the definition of 
the Riesz fractional derivative, see e.g. [36| 061 031 [331 122] for further details. This is a linear 
operator V" acting as follows 

Kf(r)=i a / -—\k\ a ^*f{k) 1 (2.7) 

where the right-hand side as to be understood in a distribution sense and f(k) denotes the 
Fourier transform. Some elementary properties are, see e.g. [22J 

V?V^ = V? +/3 (2.8) 

d 

V*=J>* = A r (2-9) 

[V> 4 ] =a0iVr 2 (2-10) 
[V° r 2 } = 2a(r ■ P )V? _2 + a(d + a - 2)V%~ 2 (2.11) 
V£./(/*r) = W a V a r f{nr) (2.12) 

We emphasise that the Riesz fractional derivative acts in many ways as a 'square root' of the 
Laplacian, see eq. (12.91) . 

Now consider the generators 

y« — -ft. 

1 -1/2 ■— U t 

Y% 2 ■= -td^ z - 2 - fir, 

Xo := - Z -td t - l -{vd r )- X - 

X 1 := (~t 2 d t -t(r-d r )-(x + Z)t)vr 2 -^r 2 
M := -/I 

Rij := r.Oj - i-jiij -/.'„ (2.13) 
Using the properties (I2.10|l2.1ip the commutators (II. ip are seen to hold true, except for 

[X 1 ^ ) }= Z —^t 2 d l V z - A S 1 (2.14) 



where 



zl id t - v; + ^ L + e - d - ± ^) (2.15) 



Hence, the generators ( 12. 13ft close into a Lie algebra on the quotient with respect to solutions 
of the equation 

S$(t, r) = (zfid t - V* + ^ (x + £ - d + l~ 2 ) ) ®(t, r) = 0. (2.16) 
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As before, the representation (I2.13P generates dynamical symmetries of the equation (12.161) 
which can be seen from the commutators 

[S, y W] = [S, W°] = [S, M ] = [5, i^] = 

2 2 

[5,X ] = [S,X 1 ] = -ztV z - 2 S. (2.17) 

Some comments are in order: 

1. because of (12. 9p . the generators ( 12. ip for an even dynamical exponent z = 2n are included 
as well, with the same invariant differential equation (I2.16p . 

2. only if one adopts the correspondence d r t— > V r between the partial derivative with 
respect to the ID spatial coordinate r and the Riesz fractional derivative V r , one has a 
relationship between the ID representation (1 1.7ft and the representation ( 12.131) for d > 1 
dimensions. 

3. the non-locality only enters into the Galilei- and special transformations, generated by 
Yf } and Xi. 

2 

4. certainly, the choice of the Riesz fractional derivative is merely motivated by its convenient 
algebraic properties. It is still unknown whether this is an appropriate choice for the 
treatment of physical systems. 

Summarising, the representation of age(d) proposed here explicitly uses generators acting non- 
locally on space. In Fourier space, the form proposed here will lead to local, but non-analytic 
generators. The special case of an even-valued dynamical exponent appears to have rather 
special and possibly non-generic properties. 

In addition, we have shown here that certain equations of motion can be said to posses the 
kind of non-local symmetry considered here, at the price of somewhat weakening the meaning 
of the term "dynamical symmetry" with respect to how it had been used for the cases where 
z = 1 or z = 2 admit generators acting locally. It remains an open problem how to reformulate 
this kind of non-local symmetry in terms of an invariance of a field-theoretic action. 

3 Co-variant two-point functions 

The two-point function 

F{t u t 2 , r u r 2 ) = rO^ia, r 2 )> (3.1) 

can be derived from the covariance under representation ( 12.131) discussed in the previous section. 
Here 

T\ = (rn, r u , -.-rid), r 2 = (r 2 i, r 2i , ...r 2d ). (3.2) 

Later, we shall mainly concentrate on the case z — 4, relevant for the physical applications we 
have in mind (see section 4), but for the moment z will be left arbitrary. As usual in conformal 
or local scale-invariance, one forms from (I2.13P two-particle operators and co-variance then 
gives the system Xi(2)F = 0. 
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1. The co- variance under space-translations requires: 

YH : (d ru + d r2i )F = (d ri - d n )F = 0, r t = r u - r 2l . (3.3) 

2 

Letting r = (n, ...r^), it follows that F = F(t±, t 2 ,r). 

2. The covariance under a "mass" transformations gives the Bargman super- select ion rule 
A*i + A*2 = and consequently only fields with opposite "masses" will give a non-trivial result 
for the two-points function. It follows that 

$i : $2 : (-//;z 2> 6), ( 3 - 4 ) 

where £i,£2 and £i,£ 2 are the scaling dimensions which characterise the corresponding non- 
stationary quasi-primary scaling operators. Consequently, the scaling operator <3>! with a posi- 
tive mass n > will be associated with the order parameter whereas the scaling operator $ 2 
with a negative mass —fi < is interpreted as the corresponding response operator <fi in the 
context of Janssen-de Dominic theory, see section 1. 



3. Next we rewrite the two-point function as 

F = F(r,v,r) , r:=t 1 -t 2 , v:=t 1 /t 2 , 

and obtain from the three co- variance conditions Y^F = for i — 1, 

2 

Xl-F = the equations (with an implicit sum over repeated indices j = 1, . . . , d) 



(3.5) 

d, X F = and 



2 T v~^l r - 2 TVdv 



TV 



-rjd r . - + 



z 

TV 



[-Tdr^-H^F 
1 „ X\+ x 2 



Td T Tjdr- — 

2 3 3 



{X2 + 6 



72-2 



2 

2 r V 



(3.6) 
(3.7) 

(3.8) 



Acting with 2 on (13. TJ) , the difficult eq. (I3.8P can be simplified to 

1 



72-2 



H 

2 -u — 



(xi - x 2 + 2ft - 2n + 2) 



x 2 - X! + 26 - z + 2) F = (3.9) 



In analogy with the one- dimensional case [23], it is clear that each of the equations ( I3.6|3.7|3.9|) 
will fix the dependence of F — F(t, v, r) on one of its variables. In fact, since the variable v 
does not enter explicitly into eqs. ( 13. 6113. 7p . these two alone will completely fix the dependence 
of F(r,v,r) on r and r such that the dependence on v factorises. Then the (z — 2)-fold 
(fractional) derivatives in (13. 9p can be dropped. We also checked that the closure condition of 
our representation is automatically satisfied, as it should be, and as seen for d = 1 before [23] . 
Isolating the usual power-law pre-factor which only depends on the single time t 2 , we find the 
following reduced form 



~(x 1 +x 2 )/z 



(v-iy 



:[(si+x 2 )/2+£i+f2-*+2] 



.l{x 2 -x 1 +2&-z+2] j ^ ^ _ 



(3.10) 
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where the last function f(r, r) has to be determined from eqs. (I3.6|l3.7p . 

4. The co- variance condition under each rotation Rij gives, for scalar quasi-primary operators 
$1.2 

(r t d rj -r,d n )F = (3.11) 
and leads to / = /(r, u), with u = r 2 as expected. 

Combining this with dynamical scaling eq. (13. 7ft . one has simply 

/ = / (ti'/V 1 ) (3.12) 

and the form of the last remaining function / has to be found from Galilei-covariance, as given 
by (E7j). 

This is the technically most demanding part of the calculation. For the sake of computational 
simplicity, we shall distinguish the cases of an even dynamical exponent z = 2n (in particular, 
we shall describe below an explicit application of the case z — 4) before we consider the case 
of generic values of z. 



3.1 Even dynamical exponent z = 2n 

Using the notation p := it^r -1 , Galilei-covariance gives, for each % — 1, . . . , d 

{rdrAr 1 + m)f{r, r) = n {2Td u A n u ~ 1 + A*) f(r, «) = U ((2n)> 2 ^ ^A^ 1 + /x) f(p) = 0. 

(3.13) 

Here we use the following notations for the various forms of the Laplacian 

^ d 2 

A r = ^u = 2dd u + Aud 2 u 

i=l * 

= ^(^+2("-D)P^^ + 2«P^)=:-^A, (3.14) 

In contrast to A" -1 which is a formal power series of A r , the explicit forms of A" -1 and A" -1 
must be calculated for any specific value of n. For example 

A 2 u g(u) = A u (A u g(u)) = (4d(d + 2)d 2 u + 16(d + 2)udl + 16w 2 ^) g(u). (3.15) 

For illustration, we solve the equation ( 13 . 1 3[) for n = 2, that is z = 4. In this case, one has 
((d + 2)d p + 2{d + 8)pd 2 p + 8p 2 c)J + |) /(p) = 0. (3.16) 
There are three independent solutions which are sought as Frobenius series 

oo 

f(p) = J2 a xP k+a > «o^O. (3.17) 

fc=0 

where the exponent a is found from the characteristic equation a (8a 2 + 2a (d — 4) + 2 — d) = 
with solutions cr = 0, |, | — |. Standard methods give a recursion relation for the coefficients 
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a m and lead to, for n = 2 or equivalently z = 4 

+A ^^(l-|,|-J-g). (3.18) 

where 0-F2 is a generalised hyper-geometric function and /o,i,2 are normalisation constants. One 
could use the explicitly known asymptotic forms of 0-^2(05 b; x) for x — > ±00 in order to identify 
the sub-space of solutions which do not diverge for \jip\ 3> 1 [561 [23], but we shall not need this 
explicitly later. 



3.2 Generic dynamical exponent 



The covariance under the generalised Galilei transformations of the representation (1 2 . 1 3 [) re- 
quires the d conditions 



(rd rj V z ~ 2 + fir J )f(r,r) = , j = l,...,d. 
This can be solved in Fourier space where 

1 



f(r,r) 



dk e lfe - r /(r, fc) 



(3.19) 



(3.20) 



which leads to the correspondences dj and V£ ^ (i|fc|) a - It follows that 
equation (13.191) takes the following form in Fourier space: 

(fiB kj +i z - 2 rk j \k\ z - 2 )f(r,k) = 0, 3 = l,...,d 
and has the the solution 



f(T,k>) = / (r)exp 

This is rewritten in the direct space as follows 
/o(r 



i 2 ~ 2 r 



Z jj, 



/(r,r) 



(27T) 



dfc exp 



i*-2 r 
ifc ■ r fe 

Z jJL 



Ml) 

(27r) d 



P -2 T 



with the following abbreviation 



Ip{r) := I dk exp [ife • 1 — /?|fe| 



(3.21) 
(3.22) 

(3.23) 
(3.24) 



In what follows, we shall admit /3 G C, but suitably restricted such that integral converges (i.e. 
Re P > would be sufficient). To evaluate the integral f)3.24p formally, introduce d- dimensional 
in spherical coordinates 



h 

kd-i 
k d 



k cos 9i 

k sin 6\ cos 6> 2 



k sin 6\ sin # 2 ... sin(^_ 2 ) cos(^_i) 
/csin #1 sin 6> 2 ... sin(0 d _ 2 ) sin(0 d _i), 



(3.25) 
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where k — \k\, 6j € [0, %) for j — 1, . . . d — 2 and 0^-1 € [0, 27r). The jacobian is 

gin d-2 &i gin d-3 ^ _ _ _ S in(0 d _ 2 ) (3.26) 



d{k,e u ...,e d - 1 ) 

Then the integral (j3.24p becomes 



pOO f"K 

Ip(r) = / dkk d - l e~ pkZ / d^e ifcrcosei sin^ 2 ^ 
Jo Jo 

x /"dfla sin^- 1 )" 2 9 2 ... r&6 d - 2 sin^- 1 )-^- 2 ) fl d „ 2 / "dfl^ 
Jo Jo Jo 



J*rcoB*i gin <i-2 



pOO p 71 

S d -i / dkk d ~ l e- pkz / d^e 1 
Jo Jo 

poo 

Jo 

pOO 

( 27T ) d / 2 r- d duu d/2 J d/2 -i(u) e^ r ~ z)uZ (3.27) 



where we used [TBI eq. (8.411.7)] in the third line in order to express the last angular integral in 
terms of a Bessel function J v {w) and S d = 2ir d / 2 /T(d/2) denotes the surface of the d- dimensional 
sphere. 

In order to display explicitly the space-time scaling of the scaling function f(r,r), we re- 
express (3 = ar such that a = i z ~ 2 / (z/i) and then have 

'<'■'> = <§^ A £) < 3 ' 28 > 

where the scaling function A is given by the above integral. Applying scale-invariance, gives an 
obvious differential equation for the last remaining function, hence /o(t) = f(o)^ d ^ z - Finally, 
use the expansion of the Bessel function J„(w) = J2^=o ^~n\r(u+n+i) ' * n or d er ^° convert the 
scaling function f(r, r) = (2n)~ d / 2 f(o)( rZ /t) 6 ^ A(arr~ z ) into a series. This gives the final result 

f(rr) f r ^ 2 )Vf 1Y> T ^ ( ^ V T3 2Q1 

f (T ' r) = /o ° vWz) n!T(n+ d ,) ) ( ] 



and where foo is a normalisation constant. Clearly, this series has an infinite radius of conver- 
gence for z > 2. 

The sought form of the two-point function F is given by eqs. (I3.10f3.29p . The special case 
of even values of z treated above is compatible with this formjj 



3 For example, in the case z = 4, one may separate the sum into even and odd terms. Repeated application 
of the duplication formula T(2k) = n-^2 2k - 1 T(k)r(k + l/2) shows that the sums arising are exactly of the kind 
found in (|3.18|) . We have also seen that the form (|3.29[) has a Fourier transformation, which is not obvious for 
all solutions in (|3.18[) . 
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4 Spherical model and field-theoretical description 



4.1 Spherical model with conserved order parameter 

The spherical model was conceived in 1953 by Berlin and Kac as an exactly solvable mathe- 
matical model for strongly interacting spins and has proved ever since to be an useful model 
for trying out more general ideas. It is usually defined in terms of real spin variable S(t, x) 
attached to each site x of the hyper-cubic lattice A C Z d and depending on time t, subject to 
the mean spherical constraint 

(^ft,*) 2 )^, (4.1) 

where N is the number of sites. The spin Hamiltonian is W = — J2( x y ) S X S V where the sum is 
over pairs of nearest neighbours. At equilibrium, this leads to a second-order phase transition 
with a critical temperature T c > for all spatial dimensions d > 2 and where the critical 
exponents have non-mean-field values for d < 4, see e.g. [28] and references therein. We 
consider here a kinetics given by a Langevin equation with a conserved order parameter (model 
B in the terminology of [25] ) 

d t S(t, x) = -Vl[V 2 x S(t, x) + }(t)S(t, x) + h(t, x)} + v (t, x), (4.2) 

where $(t) is the Lagrange multiplier fixed by the mean spherical constraint and the coupling to 
the heat bath with the critical temperature T c is described by a Gaussian noise rj of vanishing 
average and variance 

(r]{t, x)r]{t', x')) = -2T c V 2 J(t - t')5(x - x') (4.3) 

and h(t, x) is a small external magnetic field (required for the computation of the response) 
|30[ EU H3 H]. This constitutes a qualitatively reasonable model for the kinetics of phase- 
separation, for instance in alloys. 

Essentially the same equations can also be used to describe very different physical situations. 
For example, the growth of interfaces on a substrate with a conservation of particles along the 
interface is described by the Mullins- Herring equation [39] (or the equivalent Wolf- Villain model 
[55] ) which in our notation is simply given by the Langevin equation (14. 2p with 3 = 0. Several 
inequivalent variants exist: one may either take the conserved noise (14. 3p or else consider the 
non-conserved centred noise with variance (77ft, x)r](t', x')) = 2T c 5(t — t')5(x — x'), see [H] for 
details. 

Here, we shall concentrate on the Mullins- Herring model with the conserved noise (14. 3 1) . 



Since correlators and response have already been studied in detail in the literature 
Wf\ H] it is enough to quote the results. For our purposes, we merely need the full space- 
time response in the conserved spherical model for d > 4 or equivalently the Mullins-Herring 
equation for any d, which reads 

R(t,a;r)= -V— — ft - 



-(d+2)/4 



2 M / 2 7T rf / 2 r(d/4) 



0^2 



2' 4' 256(t- s) 



8T( 



dT(j+ 1 



l X(- 



, 3 d 
F 2 ( -, - + - 



2' 4 2' 256ft -s) 



(4.4) 
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We want to compare this with the a£je(<i)-covariant two-point function, as obtained in the 
previous section from the non-local representation (12.11) . with z — 4. Before, we can do this, 
however, we still have to show how to relate a dynamical symmetry of a deterministic equation 
such as (12.161) with the properties of a solution of a stochastic Langevin equation ( 14.21) . 



4.2 Field-theoretical description 

In order to compare co-variant response computed from the non-local representations and 
derived in section 3, we appeal to standard methods in non-equilibrium field-theory, see e.g. 
I5"4"] for introductions. Thereby, we shall generalise previous applications of the method 
El IH El El ED] to the present case of non-local representations. Consider the Langevin 
equation 

d t <t> = -^l (-V*0 + v{t)cf>) + V (4.5) 
with the conserved centred gaussian noise 

( V (t, r) V (t', r')> = ~^V 2 r S(t - t')S(r - r>). (4.6) 

Following [531 EH H] , the Janssen-de Dominicis action is 



J((f>,<f>) = dudR 



where 



(d u - ^ 2 R (yl - 

+ ^ JdudR0(u,R)(V 2 0(u,R)) + J init ((j),0) (4.7) 

Jinit(<f>, 4>) = \J dRdR ' ?(°» R ) <^(°» R )<f>(°' R ') ?(°' R ') ( 4 - 8 ) 

with complete analogy with non-conserved case [351 H2] • 

Averages of an observable A are defined as usual by the functional integral 

(A) = J V[<j>]V[$l A[(j>\eM-J{M)- (4-9) 

Next, one decomposes the action into a deterministic part J$ and a noise part Jb 

J$,<t>\) = Jo\iA\) + Jb\<t>\ (4-10) 

with 

1 



Jo(d>, 6) = / dudR 



du - ^r^r ~ «(«)) 



(4.1i; 



Jb = Jimt + Jth] Jth = j- J dudR <j)(u, R)(V 2 (j)(u, R)). (4.12) 

The point of this decomposition is that the action j7o is invariant under the non-local represen- 
tations constructed in section 2. We call the theory with action J7q noise free or deterministic 
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and denote the average with respect to it by (...)o- Then any average of the full theory can 
then be computed by formally expending around the noise-free theory 



(A) = (Aexp(-J b )) . (4.13) 

We can focus on the deterministic theory with v = because of the following observation: 
If (p{t, r) is a solution of equation ( [^.<5| ] with 77 = then the field 

ij>(t, r) := exp Q- jf dr u(r)V^ </>(t, r) (4.14) 

fulfils the same equation ( |^.5| ) with v = (and rj = 0). Hence, it is enough to consider first the 
problem with v = and apply the inverse of the transformation (I4.14p for treating the case 
v ^ 0, in order to implement the breaking of time-translation invariance. On the other hand, 
the case v = is relevant for d > 4 in the spherical model and for any d in Mullins-Herring 
model, which are the cases we are interested in. 

The noise-free theory is given by the equation (14. 5 j) with rj = 0, v = or by the action J Q 
with v = 0. This equation coincides with the equation (12. 5 j) in section 2 with n = 2 when 
x + £ — (d + 2)/2 = 0. If the fields 0, are indeed quasi-primary, that is they transform under 
non-local representation (12. ip of the ageing algebra agt(d), it follows that 

4> = <P(n,x,o , 4> = { -^,i) (4-15) 

are characterised with opposite "mass" parameter and scaling dimensions (x, £) corre- 

spondingly, with relation 

x + £ = x + £ = (d + 2)/2. (4.16) 

The noise-free theory then is invariant under the non-local representation (12. ip of agc(d). The 
Bargman super-selection rule [2] holds in term of the "mass" parameter \x that is 

(0 1 (t 1 ,r 1 )...0 n (t n ,r n )) o = (4.17) 

unless fii + . . . fi n = 0. On the other hand because of the Gaussian structure of the action J$ 
[1] of the noise-free theory, Wick's theorem applies and gives [SZ] 

(0l(tl,Ti)...02n(*2n,^2n))o = ^2 

All possible pairings P of 1, 2, ...2n 

(4> Pl (t Pl , r Pl )(f)p 2 (tp 2 , r P2 )) ...((f)p 2n _ 1 (tp 2n _ 1 , rp^^p^tp^ r ft J) . (4.18) 

As in the non-conserved case |l2j H], it follows that the response function is independent of 
noise, while in noise-free theory it is given by 

R(t, s-x-y) := (0(t, x)V 2 J(s, y)) = V 2 r F^(t, s; r) (4.19) 

where F^(t, s; r) is the two-point function computed in the previous section in eqs. (I3.10|3.18p ; 
after taking time-translation invariance into account, which we achieve under important as- 
sumption x + x = d, in addition to the relations (I4.16p . It follows 



R(t, s; r) = (t- sY^Vlf (^-^ = (t - sY^\f{ P ) 
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4(t - s)~ {d+2)l \(d + 2)ps0 p + 4p§a 
11 d up 



fo 0F2 



2' 2 



-(rf+2)/4 



/l 0*5 





HP 




64 


d 


3 


4 


' 2 






( 


2' 



f 1/2 Z? / 3 d , 1 ^ 



4' 64 
Id HP 



fy /2 F 2 



3 d 
2' 4 



64 



+ / 2 P oF2 (2~4' 2 ~4'"64 



(4.20) 



While the response function must be regular for r = and goes to when r — > 00 the third 
term is eliminated and the constants /q and /{ must be related. Using the known long-term 
behaviour of the hyper-geometric function [5T)| I2"3] we recover exactly the result (I4.4p . 

Therefore, we have seen that the exact result (I4.4p . rather than being interpreted as in j4] 
a consequence of a local symmetry, the closure of which as a Lie algebra is somewhat artificial 
and ad hoc, may also be seen as evidence of a non-local representation ( 12. ID of the well-known 
Lie algebra oge(d), where the closure is naturally provided by the same equation of motion 
which is under study ! 



5 Conclusions 

In trying to construct a closed set of Lie algebra generators for generalised scale-transformations 
with an arbitrary dynamical exponent z, we have been led to consider non-local representations 
of the the ageing algebra agt(d). Our construction was guided by the previously known results 
for d — 1 space dimensions [23] to generic d > 1. It turns out that the result depends on the 
value of the dynamical exponent z. 

Conceptually, we have admitted a slight extension in what is called a dynamical symme- 
try, with respect to the usual definition of that notion in the literature, see for example |41j . 
Conventionally, the infinitesimal generator X of a dynamical symmetry of the equation of mo- 
tion Sip = must satisfy [S, X] = \x<S as an operator, where is a scalar. The definition 
used here only requires this when the quotient of the space of functions ip with respect to 
the space of solutions of Sip = is taken. Such a quotient space construction appears to be 
necessary whenever the dynamical exponent z ^ 1,2. At least, this is a more elegant and 
hopefully mathematically more sober way to achieve the closure of the various generators into 
a Lie algebra, and appears to us preferable to the more traditional approach where several ad 
hoc super-selection rules had to be imposed in order to cut off the infinite set of non-closing 
operators, see [T71 141 12T1 122] . 

Several details depend on the value of z 

1. For z even, the algebra (12. ip contains the usual local generators of space-translations, 
spatial rotations, dynamical scaling and global phase shift. In addition, the further d+l 
generators of generalised Galilei-transformation and special transformations are non-local 
and are constructed with linear differential operators of order z — 1. By analogy with the 
ID case [23], we suspect that these might be interpreted as generating transformation of 
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distribution functions of the positions, rather than bona fide coordinate transformations. 
The constructed representation acts naturally as a dynamical symmetry on the solution 
space of the simple linear Schrodinger equation we considered, and should be viewed 
essentially as a toy model. The example studied here might be the first step towards 
an understanding how to use such non-local transformations in applications to the non- 
equilibrium physics of strongly interacting particles. 

We recall that in the context of interface growth with conserved dynamics, exactly the 
kind of non-local generalised Galilei-transformation we have studied here has already 
been introduced in analysing the stochastic equation (related to molecular beam epitaxy 
(mbe)) 



where r/ is a centred and conserved gaussian noise and v, A are constants [52]. Indeed, they 

show that Galilei-invariance leads to a non-trivial hyperscaling relation, expected to be 
exact and checked in the lowest orders of the loop expansion. In particular, they obtain 
z = 4 in d = 2 space dimensions [52]. Similar results exist for slightly different non-linear 
conserved growth equations [15]. We hope to return to a symmetry analysis of these 
non-linear equations in the future. In any case, the available evidence that generalised 
Galilei-invariance could survive the loop expansion is very encouraging. 

Since the non-local representations only close over solution of the same equation, we deal 
with an on shell representation]^ It would be interesting to see whether a systematic 
relationship with conditional symmetries in partial differential equations can be found, 
see [7] and references therein. 

2. For generic values of z, the generalisation from the one-dimensional case requires the 
explicit introduction of some kind of fractional derivative. For our purposes, the Riesz 
fractional derivative turned out to have the required algebraic properties. In addition, 
the result derived for the co-variant two-point function is compatible with the directly 
treatable case when z is even. We are not aware of confirmed physical applications with 
generic values of z. It remains to be seen if similar symmetries of other non-local equations 
can be found. 

On the other hand, generalisation proposed here in terms of strongly non-local Riesz frac- 
tional derivatives into the generators permits at least to close the representation on the 
solution space of a single partial differential equation. Another advantage of this realisa- 
tion is that it naturally includes the case z even, when the fractional derivatives become 
the powers of the Laplacian. It formally can be interpreted as dynamical symmetry 
algebra and may have unexpected applications which we hope to find in the future. 

Clearly, the extension of the present method to non-linear equations remains another 
open question. 
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4 Its sub-algebra, however which includes the dynamical scaling and Galilei algebra is off shell that is, it is 
closed without any conditions. 



cW> = -V 2 z/V 2 + -(V0) 2 +T) 
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